Prodigal: prokaryotic gene recognition and translation initiation site identification

<p>Abstract</p> <p>Background</p> <p>The quality of automated gene prediction in microbial organisms has improved steadily over the past decade, but there is still room for improvement. Increasing the number of correct identifications, both of genes and of the translation initiation sites for each gene, and reducing the overall number of false positives, are all desirable goals.</p> <p>Results</p> <p>With our years of experience in manually curating genomes for the Joint Genome Institute, we developed a new gene prediction algorithm called Prodigal (PROkaryotic DYnamic programming Gene-finding ALgorithm). With Prodigal, we focused specifically on the three goals of improved gene structure prediction, improved translation initiation site recognition, and reduced false positives. We compared the results of Prodigal to existing gene-finding methods to demonstrate that it met each of these objectives.</p> <p>Conclusion</p> <p>We built a fast, lightweight, open source gene prediction program called Prodigal <url>http://compbio.ornl.gov/prodigal/</url>. Prodigal achieved good results compared to existing methods, and we believe it will be a valuable asset to automated microbial annotation pipelines.</p

Assessing chemical mechanisms underlying the effects of sunflower pollen on a gut pathogen in bumble bees

Many pollinator species are declining due to a variety of interacting stressors including pathogens, sparking interest in understanding factors that could mitigate these outcomes. Diet can affect host-pathogen interactions by changing nutritional reserves or providing bioactive secondary chemicals. Recent work found that sunflower pollen (Helianthus annuus) dramatically reduced cell counts of the gut pathogen Crithidia bombi in bumble bee workers (Bombus impatiens), but the mechanism underlying this effect is unknown. Here we analyzed methanolic extracts of sunflower pollen by LC-MS and identified triscoumaroyl spermidines as the major secondary metabolite components, along with a flavonoid quercetin-3-O-hexoside and a quercetin-3-O-(6-O-malonyl)-hexoside. We then tested the effect of triscoumaroyl spermidine and rutin (as a proxy for quercetin glycosides) on Crithidia infection in B. impatiens, compared to buckwheat pollen (Fagopyrum esculentum) as a negative control and sunflower pollen as a positive control. In addition, we tested the effect of nine fatty acids from sunflower pollen individually and in combination using similar methods. Although sunflower pollen consistently reduced Crithidia relative to control pollen, none of the compounds we tested had significant effects. In addition, diet treatments did not affect mortality, or sucrose or pollen consumption. Thus, the mechanisms underlying the medicinal effect of sunflower are still unknown; future work could use bioactivity-guided fractionation to more efficiently target compounds of interest, and explore non-chemical mechanisms. Ultimately, identifying the mechanism underlying the effect of sunflower pollen on pathogens will open up new avenues for managing bee health

Genome analysis Gene and translation initiation site prediction in metagenomic sequences

ABSTRACT Motivation: Gene prediction in metagenomic sequences remains a difficult problem. Current sequencing technologies do not achieve sufficient coverage to assemble the individual genomes in a typical sample; consequently, sequencing runs produce a large number of short sequences whose exact origin is unknown. Since these sequences are usually smaller than the average length of a gene, algorithms must make predictions based on very little data. Results: We present MetaProdigal, a metagenomic version of the gene prediction program Prodigal, that can identify genes in short, anonymous coding sequences with a high degree of accuracy. The novel value of the method consists of enhanced translation initiation site identification, ability to identify sequences that use alternate genetic codes and confidence values for each gene call. We compare the results of MetaProdigal with other methods and conclude with a discussion of future improvements

MicroRNAs Form Triplexes with Double Stranded DNA at Sequence-Specific Binding Sites; a Eukaryotic Mechanism via which microRNAs Could Directly Alter Gene Expression

<div><p>MicroRNAs are important regulators of gene expression, acting primarily by binding to sequence-specific locations on already transcribed messenger RNAs (mRNA) and typically down-regulating their stability or translation. Recent studies indicate that microRNAs may also play a role in up-regulating mRNA transcription levels, although a definitive mechanism has not been established. Double-helical DNA is capable of forming triple-helical structures through Hoogsteen and reverse Hoogsteen interactions in the major groove of the duplex, and we show physical evidence (i.e., NMR, FRET, SPR) that purine or pyrimidine-rich microRNAs of appropriate length and sequence form triple-helical structures with purine-rich sequences of duplex DNA, and identify microRNA sequences that favor triplex formation. We developed an algorithm (Trident) to search genome-wide for potential triplex-forming sites and show that several mammalian and non-mammalian genomes are enriched for strong microRNA triplex binding sites. We show that those genes containing sequences favoring microRNA triplex formation are markedly enriched (3.3 fold, p<2.2 × 10⁻¹⁶) for genes whose expression is positively correlated with expression of microRNAs targeting triplex binding sequences. This work has thus revealed a new mechanism by which microRNAs could interact with gene promoter regions to modify gene transcription.</p></div

Detection of DNA-DNA and RNA-DNA triplexes by EMSA and NMR, and molecular modeling of miRNA-duplex DNA triplex.

<p><b>(A)</b> EMSA; 5’ ROX-labeled hairpin duplex DNA (0.1 μM) was incubated for 3-hrs at 22°C in the presence (lanes 2–11) or absence (lane 1) of 2.5 μM 483-opti DNA oligo, and increasing concentration (30, 60, 150 μM) of Hoogsteen bond-optimized hsa-miR-483-5p (483-opti, lanes 3–5), hsa-miR-483-5 (483, lanes 6–8), or a scrambled RNA oligo (Scramble, lanes 9–11). Duplexes and triplexes were resolved on a 20% non-denaturing acrylamide gel, and the ROX-signal visualized. Triplex of 483-opti DNA oligo and duplex DNA is readily detected (lane 2). The 483-opti RNA oligo competes with 483-opti DNA oligo for binding to duplex DNA which is evident by increased amounts of duplex DNA and decreased amounts of triplex (compare lanes 3–5 with lane 2). Hsa-miR-483-5p (483) and scrambled RNA, because of the fewer number of favorable Hoogsteen bonds, did not compete with the 483-opti DNA oligo for binding to duplex DNA (lanes 6–7 and 9–10, respectively). <b>(B-C)</b> NMR; Two-Dimensional (2D) [¹H, ¹H] TOCSY spectra of free single stranded hairpin duplex DNA (blue contours), hairpin duplex DNA combined with hsa-miR-483-5p RNA oligo (green contours; 1:1.5 ratio), and hairpin duplex DNA with single stranded DNA oligo with the same sequence as hsa-miR-483-5p (red contours; 1:1ratio). <b>(B)</b> Thymidine cross-peaks between H6 and H7 (methyl), and <b>(C)</b> cytosine cross-peaks between H5 and H6. Single stranded RNA (hsa-miR-483-5p) or single stranded DNA with hairpin duplex DNA show similar improvement in peak the intensities, and similar chemical shift perturbations/appearance of new peaks highlighted in blue boxes, suggesting that single stranded DNA and single stranded RNA of the same sequence bind to DNA duplex in a similar manner; the major differences (peaks in red boxes) are one peak among thymidine cross-peaks, showing an intermediate change (peak disappearing) with singe stranded RNA while saturated with hairpin duplex DNA, and two new peaks among cytosine cross-peaks showing much higher intensities with single stranded DNA, indicating that the latter DNA binds to duplex DNA duplex with higher binding affinity than RNA, consistent with the results obtained by EMSA. <b>(D)</b> Molecular model of hsa-miR-483-5p-DNA triplex. (I): the model of predicted miRNA and corresponding DNA duplex sequences (16 favorable Hoogsteen pairings). All predicted Hoogsteen base pairs are well maintained after removal of positional and distance restraints(II): negative control (antisense hsa-miR-483-5p) of model with 9 favorable Hoogsteen pairings. Both RNA and DNA duplex are largely twisted and nearly all predicted Hoogsteen pairings cannot be stably maintained. Residues in favor of Hoogsteen hydrogen bond formation are shown in red while the others are shown in blue.</p

MicroRNAs form triplex structures with DNA.

<p><b>(A)</b> Duplex DNA identified by genome-wide screens of binding sites was incubated in presence or absence of a synthesized hsa-miR-483-5p with a 3’ ROX label to perform a FRET assay to detect triplex formation (illustrated in <b>3B</b>). In the absence of ROX labeled hsa-miR-483-5p (<b>3A</b>, black line) a single emission peak at 520nm is observed which, with the addition of ROX labeled hsa-miR-483-5p (<b>3A</b>, red line), is diminished and a second FRET induced emission peak at 610nm is observed. <b>(C)</b> In a complementary surface plasmon resonance (SPR) based assay (illustrated in <b>3D</b>), a 3’ biotin labeled hsa-miR-483-5p was immobilized and duplex DNA was introduced in triplicate in a 2-fold dilution series starting at 20 nM.</p

Characteristics of triplex forming microRNA.

<p>The top 1 percent of Homo sapiens triplex interactions (grades 1–4) were characterized by <b>(A-C)</b> microRNA dinucleotide frequency, <b>(D)</b> microRNA length, and <b>(E-H)</b> single nucleotide frequency, and compared to these same characteristics for all human microRNAs. The percentage of purine content was the largest discriminating factor in predicting triplex formation, with the majority of binding sites having greater than 75% purine or pyrimidine content (A). Higher GC content (B), length between 21 and 25 nucleotides (D), greater than or less than average G or C content (F and G), and lower than average U content (H) also predicted triplex formation.</p

Higher expression of microRNAs forming triplex structures with duplex DNA is more frequently associated with increased gene expression.

<p>MicroRNA and mRNA expression were measured in leukemia cells (ALL) obtained at the time of diagnosis from two cohorts of patients (St. Jude Protocols Total 15 and Total 16). Genome-wide linear correlations between microRNA expression and mRNA expression calculated to form grade 1 triplex structures were assessed in each cohort separately and then a meta-analysis was performed. <b>(A)</b> The distributions of Spearman p-values for associations with positive or <b>(B)</b> negative correlations are shown. Over-representation of small p-values for positive associations was significantly enriched as compared to negative associations.</p